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Abstract 

Results for transmission through the poly{G} DNA molecule are presented. We show that (i) 
periodically arranged sodium counter-ions in close proximity to dry DNA gives rise to a new con- 
duction channel and aperiodicity in the counter-ion sequence can lead to a significant reduction in 
conduction, (ii) modification of the rise of B-DNA induces a change in the width of the transmission 
window, and (iii) specifically designed sequences are predicted to show intrinsic resonant tunneling 

behavior. 

PACS numbers: 72.10.-d, 85.65. +h, 87.15.-v, 73.22.-f 


The growing interest in DNA for nanotechnology and molecular devices stems from the 
ability of controlled growth of nucleotide sequences and its surprising conducting properties. 
Recent transport measurements on DNA strands lying between metallic electrodes have 
produced contradictory results. Fink et al 1 observed metallic behavior in A DNA using an 
electron projection microscope set-up, Kasumov et al. 2 observed superconducting behavior, 
Porath et al. 3 observed semiconducting behavior with a well defined poly{G}-poly{C} DNA 
molecule, and de Pablo el al. 4 observed insulating behavior using a scanning force microscope 
based set up. Biochemists have established a link between oxidative damage of DNA due to 
hole transport and the occurrence of genetic mutations. It has also been shown that oxidative 
damage is sequence dependent. 5 Experiments 6 ' 7 and modeling 8,9 have shown that both short 
range coherent tunneling transport and long range incoherent / hopping transport 8 10 are 
mechanisms for hole transport. 

Charge transport in DNA is a complex phenomenon because the environment of the 
DNA plays a significant role in determining ;he energy levels of nucleotides and the 
structure. 11 Transformation of DNA structure has been correlated to the pH of the buffer. 
In this letter, we model some aspects of DNA structure and environment that influence 
charge transport. 

Model and Method: The system modeled is a single strand B DNA molecule, which 
has ten bases per turn with a distance of 3.38 A between consecutive bases. We considered 
two cases. In case A, only the nucleic acids are included with an H atom termination of 
the base. For two bases there are six degrees of freedom which lead to a very large number 
of possible conformations. We calculate the change in transmission probability through an 
infinite poly{G} DNA molecule due to change ir two of these degrees of freedom: rise (the 
distance between bases) and twist (angle between consecutive bases). The coordinates of 
the molecule are generated by applying the translations and rotations given in Ref. 12 to a 
single nucleotide. In case B, the sugar phosphate backbone and nucleic acids are explicitly 
included with the structure taken from Ref. 12. Hydration shells are absent in our model, 
and fluctuations in DNA conformation do not exist as in the case of dry DNA. We consider 
a straight forward B DNA model because the structure of DNA lying on a substrate is not 
well understood. 13 In the calculations involving tie effect of the counter-ions (cations) of the 
phosphate group on charge transport, the backbone and counter-ions are both considered. 
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The location of the counter-ions are determined by energy minimization 14 of a one phosphate 
two sugar model of the DNA backbone. 

The transmission through the molecule is computed using a Green’s function frame- 
work. 15 More specifically, the method used follows Ref. 16, which is based on considering 
only the block tridiagonal matrix representation of the Hamiltonian to reduce numerical 
complexity. The diagonal blocks of the Hamiltonian represent the Hamiltonian of the bases, 
and the off-diagonal blocks represent interaction between nearest neighbor bases. Interaction 
with bases that are not immediate neighbors along the strand are neglected. In constructing 
the block tridiagonal Hamiltonian, we follow the work of Gadzy et al . 17 who gave a prescrip- 
tion to construct the block tridiagonal Hamiltonian for an extended system based on the 
Hamiltonian obtained from calculations on dimers or trimers. Such a calculation has been 
carried out in the context of the DNA by Ye and Jiang. 18 We use the methods of ref. 18 to 
construct the block tridiagonal blocks and ref. 16 to compute the transmission and density 
of states. Our calculations on dimers, along the DNA chain to obtain the Hamiltonian / 
Fock matrix 18 are performed using the B3LYP exchange-correlation functional 19 - 20 in density 
functional calculations using Gaussian 98. 14 Calculations of long DNA chains are based on 
the Fock matrices of the dimers and further self-consistency is neglected. The basis set used 
here corresponds to 6-21 G. More diffuse or complete basis sets affect our results little. 

The transmission through the molecule is computed using the following expression: 

T(E) = tr (r L G T r R G a ) . (1) 

In this equation, G r{a) represents the retarded (advanced) Green’s function and T L ( R) the 
coupling between the left (right) lead and the device. G r is the solution of the following 

equation: 

(ES-H-Z r L -Z r R )G r = 1, (2) 

where E represents the energy of the electron, H the Hamiltonian, S the overlap matrix 
and T, L (r) the self-energy. In the frame- work of nearest neighbor approximation, Tl(r) and 
E l{r) are given by: 

{ ^L(K) = H2l9L(R)Hl2’ (3) 

r l ( r ) — 27rf/2ilm [<?i(n)] #12 > 

where g r L(R) is the retarded surface Green’s function of the left (right) lead, and Im [x] 
represents the imaginary part of x. 


3 


Results: The highest occupied molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) of a nucleotide primarily lie on the bases. Hence it is believed that 
transport in native DNA occurs primarily along :he base pairs, and the backbone does not 
play a role apart from providing mechanical stability to the structure. 21 Fig. 1 corresponds 
to the variation of transmission when the rise of a B-DNA is decreased (model A). When 
this distance is decreased from the regular value of 3.38 A, the mam effect is a broadening 
of the conduction channels and a shift of the HOMO channel. Physically, the broadening 
occurs due to an increase in the hybridization between the 7r orbitals of the guanine. An 
increase of the rise (distance between consecutive bases) causes a shift and sharpening of 
the channel because the hybridization between 1 he tt orbitals of neighboring bases is now 
diminished. Energetically, the decrease of the distance between nearest neighbor bases by 
0.28 A produces an increase in the total energy of only 0.14 eV for a two base system. 
Thus, even if this conformation is energetically higher, it may be encountered in specific 
conditions. Also, it may be possible to perform experiments where DNA lying between 
contacts can be stretched as in recent experiments with carbon nanotubes. 22 Broadening 
of the HOMO transmission window is also seen when the twist angle is decreased from its 
regular value of 36° (Fig. 2) as recently found by Di Felice et al. in Ref. 23. In the absence of 
strain, the width of the HOMO transmission char nel is 370 meV (Fig. 1), which corresponds 
to an average nearest neighbor hopping strength between bases of nearly 90 meV. Changing 
the twist angle from 36° to 30° and 0° changes the width of the HOMO transmission channel 

to 270 meV and 1.07 eV respectively (Fig. 2). 

We now consider the effect of counter-ions using a DNA strand with the backbone in- 
cluded (model B). One parameter which is not well characterized in most of the transport 
experiments on DNA is its close environment. It is only very recently that a combined the- 
oretical and experimental study has pointed out the possible implications of the backbone 
with its close environment in hole transport. 24 Indeed, a DNA molecule in a saline solution 
is surrounded by a cloud of ionic species due tc the presence of a negative charge mainly 
located on one of the oxygens of the phosphate group. Measurements of current by physicists 
involves dry DNA. In dry DNA, where are the cations and hydration shells located, and is 
a hydration shell necessary or can it be avoided all together? These questions are extremely 
difficult to answer at this initial stage of the field. The answer will probably come from an 
interplay of high resolution STM/AFM measurements and modeling, which is well beyond 
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the scope of our work. Our model study here makes simple assumptions about the location 
of the cations without a hydration shell, to first understand what happens in this idealistic 
limit. As is shown below, there are very interesting issues even in this limit. We find that 
the cations have two main influences on conduction: 

(a) Transport could occur via a channel composed of cations. This would however ne- 
cessitate a periodic arrangement of cations and is studied for a few different cations. These 
results do not depend significantly on the distance of the cation from the phosphate group 
as long as it is between 2.3 and 2.8 A (distance to the nearest oxygen) and are uniformly 
placed along the chain. 

(b) Cations affect the on-site potential of the liases, and hence alter conductance through 
the bases substantially. Specifically, non uniformity in cation type could lead to a reduction 
in conductance through the bases. We have presented results for the decrease in transmission 
when a single cation is replaced or removed. 

We discuss (a) and (b) now. We compare the influence of Hydrogen, Hydronium and 
Sodium counter-ions on transmission of electrons through a poly{G} strand. The charge 
transfer from the counter-ion to the phosphate group is found using Mulliken population 
analysis to be 0.42, 0.80 and 0.60 for the cases of H, H 3 0 and Na respectively. So the Sodium 
atom has 10 + 0.40 electrons and its electronic properties are different from a fully ionized 
Sodium atom with only ten electrons, which would be akin to Neon. In the case of H and 
H 3 0, we find that the HOMO and LUMO channels are still on the base and so conduction 
would occur primarily through the base. In the calculations with Na, we find that new 
conduction channels occur in the HOMO-LUMO gap (Potassium and Lithium also show 
this feature but is not discussed here). Fig. 3 shows the plot of transmission versus energy 
with Na counter-ions. The local density of states at the Sodium atom (not shown) further 
prove that the new channel (located at —2.6 eV ) is due to the sodium atoms. That is, the 
phosphate group acts as a template for a helical wire consisting of Sodium atoms, conduction 
through which gives rise to the new channel. We have also verified that the transmission 
obtained by replacing the combined system of phosphate group, sugar and base by OH at 
the location of the phosphate groups leads to the occurrence of these conduction channels 
due to the Sodium w r ire. The above results indicate that depending on the location of the 
Fermi energy, the counter-ions of the phosphate group could play a role in determining 
the conductance. Of course, as already mentioned, this necessitates uniformly arranged 
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counter-ions along the DNA strand. 

We now present some results for the influence of more than one ionic species in the 
backbone. Experimentally, when DNA dries, we expect that there would be variation in 
both the type and location of the counter-ion. We consider a system with a chain of Sodium 
counter-ions, where a single Sodium atom is replaced by H 3 0. The mam physical effect of 
the H 3 0 counter-ion is to shift the HOMO and LUMO of the base closest to it, relative 
to other bases with a Sodium counter-ion. The transmission probability through such a 
system is shown in Fig. 3 (dashed line). The decrease in transmission over relevant energy 
ranges near the HOMO and LUMO is significant. As a result, we expect the conductance 
of a DNA strand with fluctuation in the type of counter-ions to be significantly altered 
compared to a strand with uniformly placed counter-ions. We also find that the shift in the 
energy levels due to H 3 0 (with respect to Sodium) is around 150 meV, which is comparable 
to the hybridization energy between nearest neighbor guanines. Thus, a random sequence of 
counter-ions along the backbone would eventually tend to localize the electron wave function 

along the bases. 

The principle of complementarity between bases and the ability to design arbitrary 
sequences has generated interest in DNA based electronics. Inspired by this, we have studied 
the possibility of intrinsic resonant tunneling ir DNA strands. The underlying structure 
in resonant tunneling is a quantum well region separated by two barriers. Motivated 
by the fact that guanine has a lower ionization energy than thymine, we consider the 
following two strands which are seemingly symmetric, poly{G}TGGGGGGTpoly{G} and 
poly{G}TGGGTpoly{G}. We expect the thymines to comprise a barrier and the guanines 
in between to form a quantum well. Three and six clear transmission resonances due to 
hybridization of the guanines are seen in Fig. 4. The transmission probability through a 
system with symmetric barriers should peak at unity. Interestingly, the barriers created 
by the thymines have a small asymmetry because the overlap Hamiltonians between G 
and T in {5’GT3’} and {5’TG3’} are not the same (these structures are not geometrically 
identical due to rotation of bases along the DNA helix). However, as the asymmetry is mild, 
the transmission is close to unity. Next, we consider a poly{G}TTGGGTpoly{G} strand, 
which has a much larger asymmetry in the barriers created by the thymines. Here, we find 
that while the location of the transmission resonances has not changed significantly, the 
peak values are smaller than 0.2 (inset of Fig. 4). It is noted that such resonant tunneling 
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behavior may be observable in an experiment where the Fermi energy through the DNA 
is modified by applying a gate voltage or in experiments looking at current as a function 
of the applied bias between the source and drain terminals. As sequences such as those 
used above and transport measurements are bcth experimentally feasible, observation of 
resonant tunneling in DNA may be possible in the near future. 

We now briefly discuss two aspects that have been ignored in our paper, the location of 
the Fermi level and electrostatics. Reliably computing the Fermi level of molecules connected 
to electrodes has proved to be a challenging problem because charge transfer between the 
macroscopic contacts and the molecules should be precisely determined. For a large molecule 
such as the DNA, it might be best to use values that can be infered from experiments. Refs. 3 
and 4 find that the Fermi energy lies in the FOMO-LUMO gap, which is in agreement 
with the expectation that a large molecule would tend to be more or less charge neutral 
at equilibrium. A recent interesting experiment 25 has made further progress by showing 
that the Fermi energy in poly{G}-poly{C} and poly{A}-poly{T} molecules lie closer to the 
HOMO (p-doped) and LUMO (n-doped) channels respectively. Experiments that involve 
an electrostatically coupled gate electrode have a significant advantage because the Fermi 
energy can be tuned from the HOMO band to the LUMO band in a controlled manner. 

Our assumption of neglecting self-consistency for long strands is strictly not correct in the 
presence of large applied biases and for aperiodic systems. The potential profile in these sce- 
narios will be determined by long range charge transfer. In our calculations, charge transfer 
is accurately accounted for only in the Fock matrix of the subsystem. We have verified that 
expanding the subsystem from two to three bases does not significantly affect our results. 
Self-consistent calculations to determine the fraction of voltage drop across the DNA and 
contacts are important in understanding the mechanism of current flow at finite voltage bi- 
ases but is beyond the scope of this paper. It might be of interest to note that self-consistent 
calculations in semiconducting carbon nanotubes (comparable diameter to DNA) show that 
a larger fraction of the applied voltage drops across the nanotube rather than the contacts. 26 

In conclusion, we have used a Green’s function formalism to study the transmission 
properties of poly{G} DNA molecules. Decrease of the distance between the bases from 
its regular value for B-DNA leads to a shift and change in width of the transmission 
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window. Some ionic species such as sodium could create new conduction channels. More 
importantly, we have shown that variation of ionic species in the backbone can cause a 
significant reduction in transmission through the bases. Finally, we explore the possibility 
of resonant tunneling in engineered DNA structures and find that it may be possible to 
observe intrinsic resonant tunneling in a system composed of guanine and thymine. 
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FIG. 1: Transmission versus energy for different distances between consecutive bases of a poly{G} 
DNA molecule without backbone. The solid line corresponds to a rise of 3.38 A, the dashed line 
to 3.1 A and the dotted one to 3.5 A. The inset corresponds to the linear response conductance at 
T = 298 K. 



FIG. 2: Transmission versus energy for different twist angles between consecutive bases of a 
poly{G} DNA molecule without backbone. Solid line: D = 0°, dashed line: Q = 30°. The 
inset corresponds to the linear response conductance at T = 298 K. 


10 





-6 - 4-2 0 2 

E (eV) 


FIG. 3: Transmission versus energy for a poly{G} B-DNA surrounded by Na ions (solid line) and 
when a local ion mismatch is created by replacing one Na by a H 3 O (dashed line). 
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FIG. 4: Transmission versus energy for poly{G}TGGGTpoly{G} (dashed) and 

poly{G}TGGGGGGTpoly{G} (solid). Inset: Symmetrical poly{G}TGGGTpoly{G} (dashed) 

and asymmetrical poly{G}TTGGGTpoly{G} (solid) sequences. 
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